AMS02 results support the secondary origin of cosmic ray positrons 
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We show that the recent AMS02 positron fraction measurement is consistent with a secondary 
origin for positrons, and does not require additional primary sources such as pulsars or dark matter. 
The measured positron fraction at high energy saturates the previously predicted upper bound for 
secondary production obtained by neglecting radiative losses. This coincidence, which will be 
further tested by upcoming AMS02 data at higher energy, is a compelling indication for a secondary 
source. Within the secondary model the AMS02 data imply a cosmic ray propagation time in the 



Galaxy of < 10® yr and an average traversed interstellar matter density of ■ 
to the density of the Milky Way gaseous disk, at a rigidity of 300 GV. 
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Introduction. The AMS02 experiment announced a 
new measurement of the positron fraction (ratio of to 
total flux) in Galactic cosmic rays (CRs) i. The new 
measurement extends to high energy, E ~ 350 GeV, with 
precision significantly superseding earlier experiments 0- 
|5| . The positron fraction is found to increase with energy, 
apparently saturating at e'^ /e"^ ^ 0.15 at E ^ 200 GeV. 

A rising positron fraction stands in conflict with ex- 
pectations based on popular diffusion models, assuming 
a homogeneous diffusion coefficient and a cosmic ray halo 
scale height that is independent of cosmic ray rigidity. 
This conflict has triggered numerous analyses invoking 
hypothetical primary sources for the positrons such as 
pulsars and annihilation or decay of dark matter parti- 
cles. 

In this paper we point out that the AMS02 measure- 
ment [il is in fact consistent with the simplest possible 
estimate due to the one guaranteed source: the sec- 
ondary production of e+ by the collision of high energy 
primary CRs with ambient interstellar matter (ISM) . The 
main result of this paper is contained in Fig. [TJ There, 
AMS02 je^ data (black) at high energy is seen to co- 
incide with an upper bound for secondary production 
(green, cyan uncertainty band), previously derived in flt| 
by ignoring the radiative losses of the positrons. 

In the rest of this paper we outline the derivation 
of Fig. [U explaining why the AMS02 result provides a 
strong hint for a secondary positron source. We com- 
ment on the implications of a rising positron fraction, 
that is not in conflict with a secondary source. Assuming 
secondary production, we then highlight the constraints 
imposed by the new measurement on models of OR prop- 
agation in the Galaxy. 

AMS02 and the secondary positron flux. While 
the propagation of CRs in the galaxy is poorly un- 
derstood, the expected fluxes of secondaries, such as 
positrons, are tightly constrained by the measurement of 
other secondaries, such as boron. This results from the 
fact that (i) different relativistic particles with the same 
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FIG. 1: Positron flux upper bound vs. data. Green shows 
the upper bound prediction of Ref. .IJ, divided by the 
flux measured by FERMI @]; bin widths correspond to 
the FERMI measurement uncertainties. Cyan shows an 
estimated calculation uncertainty on the flux upper bound. 
AMS02 [2], early 3 and subsequent ^41] PAMELA data are 
shown in black, blue and magenta. Note that extending the 
calculation to _E > 300 GeV (shaded region on the right) 
involves an extrapolation of B/C data; see text. 



rigidity (momentum divided by charge) propagate in a 
magnetic fleld in the same way, regardless of the mag- 
netic field configuration; and (ii) the production rates of 
all secondaries are correlated in a calculable manner. 

The measured number densities of stable secondary 
CR nuclei are proportional to their net local production 
rate and are thus well described by 
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where <Jj^i is the decayed spallation cross section of the 
parent nucleus j into the secondary i per ISM nucleon 
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and CTi is the cross section for destruction of i per ISM 
nucleon. The fframmaee Xesc, defined by Eq. (P), param- 
eterizes the column density of target material traversed 
by the CRs and is the same for all species. The value of 
Xcsc is measured to be 
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30% in the range 



with different fits varying by 
10 GeV < E/Z < 100 GeV 0^b|. 

Eq. ([T]) docs not capture the effect of energy loss during 
propagation. This means that it cannot be directly ap- 
plied to positrons, since they are subject to Synchrotron 
and inverse- Compton (IC) losses. Nevertheless, it was 
realized in [ij that Eq. (P) provides a robust upper limit 
to the positron flux, given that radiative losses can only 
decrease the flux of the steep positron spectrum. This 
upper limit is model independent, derived from data and 
requires no free parameters. 

Positron fraction measurements by AMS02 [2] and 
PAMELA [1, 3] are compared to the upper bound of 
Eq. (H)) in Fig.[T] Since the experiments did not report the 
positron flux, but only the / fraction, we present the 
upper bound by dividing the theoretical e"*" flux by the 
total FERMI flux measurement Q . As demonstrated 
in the figure, the upper limit is not violated by the new 
AMS02 data. This means that the data is consistent with 
secondary positrons. Moreover, at high energy the mea- 
sured positron fraction saturates within the secondary 
limit, previously predicted in [l|. This coincidence, while 
yet to be further tested by future AMS02 data at higher 
energy, is a compelling hint for a secondary source. 

It is worthwhile to compare this result to models in- 
voking new primary sources such as pulsars or dark mat- 
ter. In such models, ad-hoc tuning of free parameters is 
required to account for the positron fraction saturating 
at ~ 0.15 for E > 200 GeV. The distinction between 
the secondary and primary models is even more trans- 
parent when considering the absolute e"*" fiux. In con- 
trast to the /e^ fraction, that has a limited dynami- 
cal range, the e+ fiux due to primary sources could well 
have been orders of magnitude below or above the sec- 
ondary bound. We know of no intrinsic scale, and thus 
of no reason, in any of the primary injection models sug- 
gested in the literature, for the positron flux to lie close to 
the data-driven secondary bound, throughout the range 
£; - 10 - 300 GeV. 

A p consistency check, future tests of the sec- 
ondary model, and calculation uncertainties. A 

test of the validity of our calculations is presented in 
Fig. [21 where the measured flux of secondary antipro- 
tons [111] , that are produced in the same interactions as 
secondary positrons, is compared to the flux obtained 
from Eq. ([T}. As seen in the flgure, our calculation is 
consistent with the observations. 
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FIG. 2: PAMELA p/p data [TJl vs. the secondary source 
prediction of Eq. {T]). Cyan shows an estimated calculation 
uncertainty on the secondary prediction. Note that extending 
the calculation to E > 300 GeV (shaded region on the right) 
involves an extrapolation of B/C data; see text. 



The secondary source hypothesis will be further tested 
with upcoming AMS02 measurements of the e"*" and p 
flux at higher energy, up to the TeV range [12[. A po- 
tentially useful independent check, though complicated 
by systematic uncertainties, can be done by analyzing 
the elemental ratios of nuclei having a radioactive iso- 
tope component with a rest frame lifetime of the order of 
1 Myr, including Be/B, Cl/Ar and Al/Mg at high rigid- 
ity similar to the cooling time of the positrons [l|, [iSl, [iJ] . 
A more straightforward check, limited however to E/Z < 
10 GeV, will come from directly measuring the isotopic 
ratio i°Be/9Be. 

We now comment on the systematic uncertainties in- 
volved in computing the e"*" upper bound and the p flux. 
First, we caution that extending the e"*" and p calcu- 
lations to high energies above 300 GeV, as in Figs. [1] 
and [21 involves extrapolation. Figs. [H and [2] were com- 
puted using the CR grammage of Eq. extracted 
from analyses of HE AO 3 data and tested only up to 
E/Z 300 GeV ^]. While there are some indica- 
tions that the same parametrization of Xcsc may hold 
to E/Z ^ TeV [lij, the case is not settled with hints 
of spectral hardening reported in IB, III- Note that a 
hardening in Xcac should lead to hardening of all secon- 
daries including e"*" and p. We shade in cyan the rigidity 
range in Figs. [U and [21 where Eq. ([2|) is being extrapo- 
lated beyond the current determination of Aosc- A robust 
secondary prediction in this range will be possible with 
improved boron/carbon (B/C) measurements, expected 



from AMS02 [12 1. 



Additional systematic uncertainties, estimated roughly 
by 50% for both the e+ and p calculations, are denoted 
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by the cyan bands in Figs. [T]and[2l The main potential 
sources of error are these: 

(i) Different cross section parameterizations for hadron 
production in pp and pA colhsions, differ by energy- 
dependent factors in the order of tens of percent. The 
difficulty is the inapplicability of perturbative calcula- 
tions, together with the scarcity of accelerator data for 
soft charged hadron production, particularly at high ra- 
pidity. Resolving this ambiguity is beyond the scope of 
the current paper. Here we follow the same calculation 
done in to which we refer the reader for more details. 

(ii) We expect Eq. H]) to only apply to 10% accuracy, 
which is roughly the level at which the assumption of neg- 
ligible energy change during propagation can be expected 
to hold for stable secondary nuclei. We also estimate 
about 30% uncertainty for Xesc at 100-300 GeV/nuc. 

(iii) The primary CR nuclei flux and composition at 
the 0.1-10 TeV/nuc range, responsible for ^ 10 — 100 GeV 
e"*" and p production, are still somewhat uncertain 




Existing measurements at the relevant range |19l - l21j 
differ systematically by 20-30%. In our analysis we 
adopt a proton flux interpolating the published PAMELA 
data pj] below a TeV and the CREAM data [21] up to 
E = 100 TeV. 

Finally, we comment that the combination of AMS02 
e+/e=^, FERMI e=^, and PAMELA p and He measure- 
ments, required here in order to compare the secondary 
predictions to data, is prone to systematic errors. It 
would be preferable to use direct measurements of all 
of these ingredients by AMS02. 

On a positron fraction rising with energy. Due 

to synchrotron and IC energy losses, the positron flux is 
suppressed, compared to the upper bound, by an energy 
dependent factor /e+ < 1. /g+ should increase mono- 
tonically as a function of icooi/^csc, where tcooi is the e+ 
radiative cooling time and tesc is the mean propagation 
time. In the limit tcooi/iesc ^ 1, we expect /e+ 1. 

The claims in the literature, that the increase with 
energy of the positron fraction is inconsistent with a sec- 
ondary origin, are based on two lines of reasoning, neither 
of which is supported by data (see [H for a detailed dis- 
cussion). The first line of reasoning assumes that (i) pri- 
mary p and e~ have the same production spectrum, and 
(ii) primary e~ and secondary suffer the same energy 
losses. Both (i) and (ii) are unsubstantiated. It is plau- 
sible that primary e~ suffer additional energy loss at the 
primary CR sources and that the injected e~ spectrum 
is different than that of the protons. The second line of 
reasoning adopts some specific propagation model, which 
leads to tcooi/iosc (and so to /e+) that decreases with en- 
ergy. Such behavior of icooi/^osc and fg+ can be modified 
in alternative models. 

This discussion makes clear that given the current 
AMS02 data, depicted in Fig. [H the call for primary 
sources is unconvincing. Since both tcooi and icsc are 
neither directly measured nor reliably calculable, the en- 



ergy dependence of /e+ , and hence of the positron frac- 
tion, can not be reliably predicted. The positron data 
should be regarded as a first direct measurement of fe+ , 
with interesting implications for the times scales tesc and 
^cooi (see e.g. fiil, |23|, and more recently [ll. [l3]). 
Interpretation: constraints on CR propagation. 
In the rest of this paper we assume that the positron 
flux is of secondary origin, and proceed to deduce new 
constraints on CR propagation. 

The secondary model allows us to quantify the amount 
by which the positron flux is suppressed by propagation 
energy loss, based on the observations. The suppression 
factor /g+ is given by the ratio between the observed e~^ 
flux to the calculated upper bound. This corresponds, in 
Fig. [H to the ratio of the black data to the green curve. 
We now analyze the constraints arising from Fig. [1] 

1. CR propagation time. If we ignore Klein- Nishina 
corrections (see discussion below), then Fig. [T] implies 
that: 
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Here Ut is the time-averaged total electromagnetic en- 
ergy density in the propagation region. Note that it is 
natural to expect that Ut should depend on CR rigidity. 
Thus Ut should be understood as function of E, though 
we omit the explicit dependence for clarity of notation. 

One irreducible source for energy dependence in the ef- 
fective value of Ut comes from Klein-Nishina corrections, 
that are neglected in Eqs. ^ and (jlj which were derived 
in the Thomson limit. The Thomson limit is not a good 
approximation for 20-200 GeV positrons if Ut contains 
a significant UV component, as appears to be the case 
from estimates of the Milky Way radiation field [23|. In 
that case, the effective radiation energy density for an 
^ 200 GeV e"*" can be significantly lower than that for an 
^ 10 GeV one (see e.g. [25|). In Fig. |3] we plot the cool- 
ing time tcooi for electrons and positrons under different 
assumptions for Ut- Smooth lines set the UV compo- 
nent to zero. Dashed lines show varying amounts of UV 
light having a black body spectrum with a temperature 
T = 6000 K. 

2. The mean ISM density of the CR halo. We can 
now estimate the mean ISM density traversed by CRs. 
Using Eq. © together with Eqs. (O and (g]), we find 
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FIG. 3: Cooling time tc for e radiative losses, as function 
of energy, for different assumptions regarding the total 
electromagnetic energy density and its UV component. 



assuming an ISM composition of 90%H+10%He by num- 
ber. 

Eqs. (O and ^ suggest that the confinement volume 
of CRs decreases with increasing CR rigidity, to the ex- 
tent that CRs a.t E/Z ^ 200 GeV spend much of their 
propagation time within the thin Galactic HI disc, with 
a scale height h ~ 200 pc, while CRs at E/Z ~ 10 GeV 
probe a larger halo. These are not robust conclusions, 
however. For example, if a significant fraction of the 
grammage Xcsc is accumulated during a short time in 
relatively dense regions, e.g. near the CR source, then 
the halo could be larger. Significant energy dependence 
in Ut could further affect the interpretation. For ex- 
ample, Ut oc E~'~'-^ (inspired by the CR grammage 
Xcsc oc E~°'^) would allow for a rigidity independent 
niSM- 

Conclusions. The positron fraction measured recently 
by AMS02 is consistent with the upper bound predicted 
in Ref. 1], assuming a secondary source. Upcoming 
AMS02 measurements of the e+ and p fiux at yet higher 
energies will continue to test the model. 

At the highest measurement energy, the positron fiux 
saturates the upper bound, and throughout the measure- 
ment range it is never smaller than a factor of /e+ ^ 0.3 
compared to it. We find this to be a compelling hint for 
a secondary source. Considering hypothetical primary 
sources such as pulsars or dark matter, we know of no 
intrinsic scale in these models that would fix the positron 
fiux at this particular range. 

Interpreted under the secondary source hypothesis, 
the positron data places interesting constraints on the 
propagation time of CRs ai E/Z ^ 10 - 300 GeV, 



that we roughly summarize by tesc(10 GeV) > 30 Myr 
and tesc(200 GeV) < Myr. The constraint on igsc at 
E/Z > 100 GeV is obtained by the new positron data, 
with no direct counterpart in earlier CR data. The con- 
straint at E/Z ~ 10 GeV is consistent, within uncer- 
tainties, with measurements of the elemental ratios of 
radioisotopes [H, [13, 14 1. 

Using the measured CR grammage together with the 
new constraints on tcac, we derive the mean ISM par- 
ticle density in the propagation region of high energy 
CRs, nisM > 0.6 cm-3 for E/Z = 200 GeV. This re- 
sult for fijsM is comparable to the mean ISM density in 
the Milky Way HI disc. At E/Z ^ 10 GeV we find a 
smaller mean density, fiisM < 0.15 cm~'^. Put together, 
these numbers could mean that the scale height of the 
CR halo decreases with increasing CR rigidity (however, 
see discussion following Eqs. [5] and [5] for caveats in this 
interpretation). 
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